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Cy I We study the influence of low levels of metal enrichment on the cooling 

and collapse of ionized gas in small protogalactic halos using three-dimensional, 

OO I smoothed particle hydrodynamics simulations. Our initial conditions represent 

^ ' protogalaxies forming within a fossil H ll region, a previously ionized H ll region 

^ ' which has not yet had time to cool and recombine. Prior to cosmological reioniza- 

^ ' tion, such regions should be relatively common, since the characteristic lifetimes 



of the likely ionizing sources are significantly shorter than a Hubble time. We 
show that in these regions, H2 is the dominant and most effective coolant, and 
that it is the amount of H2 formed that determines whether or not the gas can 



O 

in 

O ■ collapse and form stars. 

^ ■ At the low metallicities {Z < 10^^ Zq) thought to be associated with the 

I ' transition from Population III to early Population II star formation, metal-line 

^ • cooling has an almost negligible effect on the evolution of low-density gas, altering 

^ ■ the density and temperature evolution of the gas by less than 1% compared to 

the metal-free case at densities below 1 cm~^ and temperatures above 2000 K. 
■ Although there is evidence that metal-line cooling becomes more effective at 

' higher density, we find no significant differences in behaviour from the metal-free 

^ case at any density below our sink particle creation threshold at n = 500 cm~^. 

Increasing the metallicity also increases the importance of metal-line cooling, 
but it does not significantly affect the dynamical evolution of the low-density gas 
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until Z 0.1 Zq. This result holds regardless of whether or not an ultraviolet 
background is present. 

Subject headings: galaxies: formation — molecular processes — stars: formation 



1. Introduction 



Within the framework of cold dark matter (CDM) cosmology the formation of structure 
proceeds in a hierarchical fashion. At high redshifts, low-mass halos with virial temper- 
atures less than ~ 10^ K are abundant. The cooling of primordial gas in these halos is 
regulated by molecular hydrogen, as H2 is the only coolant present in significant quantities 
that r emains effective at te m peratures below 10'^ K (ISaslaw fc ZipoyI 119671 : [Peebles fc Dicke 
19681 : iMatsuda et al.lll969l ). iTegmark et al.l (119971 ) developed analytic methods to model 



early baryonic collapse via H2 cooling. Numerical studies of the formation of primordial 
gas clouds and the first stars indicate that this process likely began as early as z ~ 30 



(Abel et al. 


2002; 


Bromm et al. 


2002 


)• 



develop a semi-analytic model based on the Tegmark et al. (1997) methods and included 
the effects of dynamical heating caused by the thermalization of kinetic energy of infall into 
a deepening potential. Both approaches suggest that only gas in halos more massive than 
sor rie critical mass M n rit wil l cool effectively. Much o f this work h a s rece ntly been reviewed 



bv iBromm fc LarsonI (120041 ) . ICiardi fc Ferraral (l2005f ). and iGloverl (120051 ). 



Population III stars are the first potential producers of UV photons that can contribute 
to the reionization process and are the first producers of the metals required for the formation 
of population II stars. An important question is whether small protogalaxies that formed 
within the relic H ll regions left by these first stars could form stars themselves, or whether 
the elevated temperatures and fractional ionizations found in these regions s uppressed star 
forma tion until larger protogalaxies formed. In a recent analytical study, lOh fc Haiman 
(120031 ) argue that the first stars injected sufficient entropy into the early intergalactic medium 
(IGM) by photoionization heating and supernova explosions to strongly suppress further star 
formation inside low-mass protogalactic halos in th e affected reg i ons. On the other hand, 
previous numerical work by lRicotti et al.l ( 20021 ) and lO'Shea et al.l ( 20051 ) . as well as our own 
simulations of hot, ionized gas in small protogalactic halos, show that gravitational collapse 
mediated by H2 cooling remains possible in such regions, within protogalaxies with masses 
above a certain threshold, implying that in these systems star formation is not strongly 
suppressed. 



A further form of negative feedback that must be taken into account comes from UV 
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photons within the Lyman- Werner bands of H2, which are produced in large numbers by 
massive population III stars and which can photodissociate H2, thereby quenc hing molecular 



hydrogen cooling and delay i ng the cooling and collapse of the primordia l gas (IHaiman et al. 



2000 



Machacek et al.ll200ll : iGlover fc Brand! 1200 ll : I Johnson et all 120071 ) 



Metals produced by the first stars will also be injected into some fraction of the ionized 
volume, and so will be present in gas falling into new or existing protogalactic halos. The 
question then arises as to how this low level of metal enrichment affects the ability of the 
gas to cool and collapse. iBromm et al.l (I2OOII ) studied the collapse of cold, metal-enriched 
gas and argued that there exists a critical metallicity Zcrit below which the gas fails to 
undergo continued collapse and fragmentation. However, their simulations did not allow 
for H2 formation and H2 cooling, which is expected to play an important role, unless the 
external UV background is extremely strong. Also, by starting with cold gas, they implicitly 
assume that no extra entropy or energy has been added to the gas during its enrichment; 
although, as lOh fc HaimanI (120031 ) have shown, this is unlikely to be the case. We reexamine 
this question using more appropriate initial conditions and a more detailed treatment of the 
cooling and chemistry of the metal-poor gas. 

We present simulations showing that metals play a dominant role in regulating the 
cooling of ionized gas in small protogalactic halos only for metallicities Z ~ 0.1 or 
greater. At lower metallicities, the metals have little effect and the question of whether or 
not gas in a particular halo collapses is entirely determined by the amount of H2 forming in 
that halo. We reach the same conclusion regardless of whether or not the gas is illuminated 
by a Lyman- Werner UV background, as in the absence of H2 the metals themselves do not 
provide enough cooling to allow the gas to collapse. 



2. Simulations 



2.1. Numerical Method 



To help us to assess the influence of metals on the cooling and collapse of gas in small 
protogalactic halos, we performed a number of numerical simulations. During collapse, gas 
increases in density by several orders of magnitude, and so is best simulated by a numerical 
method with a high dynamical range. We therefore chose smoothed particle hydrodynamics 
(SPH). Excellent overviews of the method, its numerical implementation, and some of its 
applications are given in reviews by iBeng (11990) and Monaghan(1992 , 2005). We use version 
1 of the parallel code GADGET, designed by ISpringel et al.l (I2OOII ). SPH is a Lagrangian 
method for simulating astrophysical flows, in which the fluid is represented by an ensemble of 
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particles, with flow quantities at a particular point obtained by averaging over an appropriate 
subset of neighboring SPH particles. As particle time steps in an SPH code are constrained 
by the Courant condition, they grow increasingly short as more particles cluster in high- 
density regions. Replacing dense cores with artificial sink particles can therefore lead to 
a considerable increase in computational performance, allowing the dynamical evolution of 
the gas to be followed over many free-fall tir nes. In the runs p resented here we introduce 
sink particles according to the prescription of iBate et al.l (Il995l ). Sink particles are created 
once the density rises above 500 cm~^ and are endowed with an accretion radius of 5 pc. On 
every time step, any gas within this accretion radius that is gravitationally bound to the 
sink particle is accreted by it. The design and imp lementation of our sink particle algorithm 
is discussed in more detail in lJappsen et al.l (120051 ). 



2.2. Chemistry and Cooling 

We have further modified GADGET to allow us to follow the nonequilibrium chemistry 
of the major coolants in both primordial and low-metallicity gas. Our chemical model is 
a simplifi ed version of the genera l model for low-density, low-metallicity gas presented in 



Paper I (iGlover &: JappsenI 120071 ). In our simplified network, we follow the chemistry of 
twelve separate species: free electrons, H^, H^, H, H^, II2, C, C"*", O, O"*", Si, and Si^. The 
more detailed model in Paper I also follows the chemistry of He, He+, Si"*""^, D, D"*", and HD. 
However, the helium and Si"*"^ chemistry are important primarily at temperatures T>10^K, 
or if the gas is partially ionized by X-rays or cosmic rays, while the deuterium chemistry is 
only significant at temperatures T < 200K, which are only reached in a few of our simulations, 
and so the neglect of these additional reactions reduces the computational cost imposed by 
the chemistry without significantly affecting our main results. The other simplification we 
made here compared to the model of Paper I is the neglect of photochemical reactions 
involving photons more energetic than 13.6 eV; specifically, we omit the photoionization of 
H, He, H2, O, and Si"*" (reactions 48, 50, 54, 57, and 59 in Paper I). We omit these reactions 
because we assume that during the simulations, our model protogalaxies are shielded by the 
neutral hydrogen in the intergalactic medium from any external sources of ionizing radiation. 
A full list of reactions included in the chemical network used in the simulations presented in 
this paper is given in Table [H Further details (including a discussion of why these particular 
reactions were chosen and full details of the rate coefficients used) are given in Paper I. 

In implementing this chemical network within GADGET, we made the decision to follow 
the nonequilibrium abundances of H2, H+, C^, O^, and Si^ and use conservation laws for 
charge and element abundance to track H, e~, C, O, and Si. The remaining two species, the 
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H~ and ions, have very short chemical equihbrium timescales and hence were assumed, 
for simphcity, to reach equihbrium instantaneously. To compute the rate of photochemical 
reactions (such as the photodissociation of H2) in runs in which an ultraviolet background 
is present, we assume a background spectrum that has the shap e of a 10^ K black body 
(as should be typical of the brightest Population III stars; see e.g. ICojazzi et al.ll2000l ). cut 
off at energies greater than 13.6 eV to account for absorption by neutral hydrogen in the 
IGM. We quantify the strength of the background by fixing the flux at the Lyman limit: 
J{^a) = 10^^^ J21 erg s~^ cm^^ Hz~^ sr^^. A list of the resulting photochemical rates is given 
in Table [21 computed using the cross sections listed in Paper I. 

If sufficient H2 forms within the protogalaxy, it will begin to self-shield, reducing the 
effective photodissociation rate. An exact treatment of the effects of self-shielding is com- 
putationally infeasible, as it would require us to solve for the full spatial, angular, and 
frequency dependence of the radiation field at every time step. Instead, we have chosen to 
incorporate it in an approximate manner. We assume that the dominant contribution to the 
self- shielding at a given point in the protogalaxy comes from gas close to that point, and so 
we only include the contribution to the self-shielding that comes from nearby H2. 

Finally, we assume that ionization from X-rays or cosmic rays is negligible. Previous 
work suggests that even if a low level of ionization from such sources is present, it will not 



have a major effect on the outcome of the collapse (iGlover &: Brandll2003l : iMachacek et al. 



2OO3I ). In any case, the initial fractional ionization in our simulations is very much higher 
than could be produced by a realistic flux of X-rays or cosmic rays. 

A second major modiflcation that we have made to the GADGET code is a treatment of 
radiative heating and cooling. In primordial gas, we include cooling from three main sources: 
electron impact excitation of atomic hydrogen (Lya cooling), which is effective only above 
about 8000 K; rotational and vibrational excitation of Ho ; and Compton cooling. Rates for 



Lya cooling and Compton cooling were t aken from ICenI ( 



cooling we used a cooling function from iLe Bourlot et al 



9921). w hile for H2 rovibrational 
(I1999I ). A number of different 



parameterizations of the H2 cooling function ha ve been used in the astr ophysical literature, 
and so we show for reference in Figure [T] how the iLe Bourlot et al.l ( 119991 ) rate compares with 
various of the other rates. 

In addition, in cases where the gas is metal enriched, we compute the amount of flne 
structure cooling coming from C, C+, O, Si, and Si+. In this calculation, we include the effect 
of excitation of the flne-structure lines by the CMB, which prevents th e gas from cooling 
below TcMB- We also include heating from H2 photodissoci ation, following iBlack fc Dalgarno 
( ll977l ). and from the ultraviolet pumping of II2, following iBurton et al. Il990h . 
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We do not include cooling from HD, as this is only significant compared to H2 at tem- 
peratures below T ~ 200 K. Dense gas which reaches these temperatures in our simulations 
is very soon thereafter accreted by a sink particle, at which point its further temperature 
evolution cannot be followed. We therefore do not expect that the lack of HD cooling in 
our models will significantly affect our conclusions. We also omit cooling from other pri- 
mordial molecules and molecular ions (e.g. LiH and Hjj"), as these are unimportant at the 
densities and temperatures encountered in our simulations. Heavy molecules such as CO 
are also omitted, as they form in significant abundances only at densities greater than those 
considered in this paper, as has already been discussed at length in Paper I. A brief list 
of the thermal processes included in our simulations is given in Table [Sj further details, in 
particular of our treatment of fine-structure cooling, can be found in Paper I. 



2.3. Initial Conditions 

2.3.1. Initial temperature and density distribution 

We initialize each of our simulations with gas that is hot (Tg = 10^ K) and fully ionized. 
The physical situatio n that these ini t ial co nditions are intended to represent is a protogalaxy 



forming within what lOh fc Haimad (120031 ) term a fossil"Hll region, an Hll region surround- 
ing an ionizing source that has switched off, containing gas that has not yet had time to 
cool and recombine. Prior to cosmological reionization, such regions should be relatively 
common, since the characteristic lifetime of the likely ionizing sources - massive Population 
III stars and/or active galactic nuclei - are significantly shorter than the Hubble time. 

The initial temperature of gas in such a fossil H ll region should, strictly speaking, depend 
on various factors such as the spectrum and strength of the ionizing source and the density 
of the gas. However, since the gas cools rapidly through Lya cooling at the start of our 
simulations, and since the temperature at which Lya cooling becomes ineffective does not 
depend on the initial temperature of the gas, we do not expect our results to be sensitive to 
small changes in the initial value of Tg. 

For our purposes, we do not need to follow the assembly history of the dark matter 
halo in which the protogalaxy resides, or to study the response of the dark matter halo 
to the cooling of the gas. Therefore, we choose to model the influence of the dark matter 
halo (hereafter simply the 'halo') by using a flxed background potential. To construct this 
potential, we assum e that the halo is spherically symmetric, with the density proflle of 



Navarro et all fll997l ): 

P"'^^''' ~ r/r,(l + r/r,)2' 
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where is a scale radius, 5c is a characteristic (dimensionless) density, and pcrit = SH'^/SnG 
is the critical density for closure. Note th at we use a value for the Hubb le constant of 
-f^o = 72kms~^ Mpc~^ ( ISpergel et al.ll2003l ). Following iNavarro et al. we calculate 



the characteristic density and the scale radius from a given redshift and dark halo mass. 
We truncate the halo at the radius at which the value of pdm given by equation [1] equals 
the cosmological background density at the beginning of the simulation. To simplify the 
discussion of our simulation results, we take as a fiducial example a halo with a total mass 
Mtot = 7.8 X 10^ Mq and an initial redshift Zi = 25. For this example halo, the virial 
temperature Tvir = 1900 K, the virial radius r^ir = 0.1 kpc, and the truncation radius rt = 
0.49kpc, where both radii are given in physical units. The scale radius of this example halo 
is 29 pc, and the full computational volume is a box of side length 1 kpc. In addition to this 
fiducial example, we also studied cooling and collapse in halos with various different masses 
in the range 4.0 X 10"^ Mq < Mtot < 1.2 X 10^ Mq, at redshifts in the range 20 < z < 30. A 
complete list of runs is given in Table HI 

We use periodic boundary conditions for the hydrodynamic part of the force calculations 
to keep the gas bound within the computational volume. The self-gravity of the gas and 
the gravitational force exerted by the dark matter potential are not calculated periodically, 
since we assume that other dark matter halos and their gas content are distant enough to 
neglect their gravitational influence. 

As we are dealing with a gravitationally bound system (the dark matter halo) that has 
decoupled from the Hubble flow, we do not include the effects of cosmological expansion, 
and we perform our simulations using physical rather than co-moving coordinates. This 
simplification should be reasonably accurate for gas within the virial radius of the halo, 
although it will break down at larger radii. However, the behaviour of the gas at radii 
r ^ r^ij. will have little effect on the evolution of the gas in the central regions of the halo, 
and since it is the latter that we are primarily concerned with in this paper, this simplification 
should be adequate for our current purposes. 

We begin our simulations with a uniform distribution of gas with an initial density Pg, 
taken to be equal to the cosmological background density, and then allow the gas to relax 
isothermally until it reaches hydrostatic equilibrium. Note that this initial phase of the 
simulation is merely a convenient way of generating the appropriate initial conditions for 
the simulation proper, and so we do not include the effects of chemical evolution or cooling 
during this phase. The mass of gas present in our simulation was taken to be a fraction 
^b/^dm of the total mass of dark matter, where the dark matter density Q^m = — ^b, 
and where fib is the baryon d ensity and flm is th e matter density. We take values for the 



cosmological parameters from ISpergel et al.l (l2003l ) of fib = 0.047 and fl^i = 0.29, giving us 
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a total gas mass of Mg = 0.19 M^m- Here, M^m is the sum of the halo mass and the mass of 
the dark matter background in the simulated volume. For our fiducial example, and for our 
example, Mdm = 1-84 x 10^ Mq. Therefore, in this case Mg = 3.5 x 10^ Mq. The number 
of SPH particles used to represent the gas depends on the mass of the halo. We selected an 
appropriate particle number and particle mass by requiring that the number of particles used, 
A^, and the particle mass, Mpart, satisfied N > 10^ and Mpart < 5.0 Mq respectively. Our 
SPH smoothing kernel encompasses approximately 40 particles and since we need twice thi s 



number in order to properly resolve gravitationally bound objects ( iBate fc BurkertI 119971 ). 
our mass resolution, in the worst case, is Mres — 80 mpart = 400 Mq. In some of the smaller 
halos we simulate, our mass resolution is significantly better. The values of and Mpart used 
for each of the different halo masses are summarized in Table O along with the corresponding 
values of Mres and also the virial temperatures of the various halos. 

At our sink creation density, Mj exceeds Mres as long as T > 30 K. This temperature is 
much smaller th a n tha t reached by the gas in our simulations, and so we always satisfy the 



Bate fc BurkertI (119971 ) resolution criterion by a comfortable margin in gas which has not 



been accreted by a sink. 



2.3.2. Metallicity 



We studied the evolution both of zero-metallicity, primordial gas and of metal-enriched 
gas. The range of metallicities studied was 10"'^ < Z/Zq < 1.0. The value used in any 
particular simulation is summarized in Table |H In particular, we performed a large number 
of simulations with Z = 10~^ Zq. This value is an upper limit d erived from Q SO absorption- 
line studies of the low column density Lya forest at 2; ~ 3 (jPettinil Il999l ). Estimates of 
the globally-averaged metallicity produced by t he sources responsible fo r reionization are 
also typically of the order of 10~^ Zq (see e.g. iRicotti &: Ostrikerl |2004| ). Moreover, the 
fact that several hundred stars with [Fe/H] < —3 have been observed in the Galactic hal o 
provides solid evidence that solar mass stars can form at this metallicity (jchristlieb 2003). 
We originally intended to examine the behaviour of the gas at even lower metallicities, but 
these additional simulations were rendered unnecessary by the fact that even at Z = 10~^ Zq 
metal cooling proved to be of no importance at the low gas densities studied in this paper. 

In our simulations of metal-enriched gas, we assume that mixing is efficient and that 
the metals are spread out uniformly throughout the computational domain. We also assume 
that the relative abundances of the various metals in the enriched gas are the same as in 
solar metallicity gas; given the wide scatter that exists in observational determinations and 
theoretical predictions of abundance ratios in very low-metallicity gas, this seems to us to 
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be the most conservative assumption. However, variations in the relative abundances of an 
order of magnitude or less will not significantly alter our results. 

As we discuss in Paper I, substantial uncertainties exist concerning the presence, abun- 
dance, size distribution, and composition of dust grains in low-metallicity gas at high redshift. 
An exhaustive examination of the various possibilities would require an unfeasibly large num- 
ber of simulations to be performed, and so we instead make the simplifying assumption that 
if dust is present, it has similar properties to Milky Way dust, but with an abundance that 
is a factor of Zj Zq smaller than the Galactic value. To allow us to investigate the relative 
importance of the dust and the gas-phase metals, we perform simulations both with and 
without dust. 



2.3.3. Other details 



We also investigate the influence of a U V background. By z = 25, a. considerabl e 



UV background may already have developed (IHaiman et al.l l200Clt iGlover fc Brandl l2003l ). 
particularly if cosmological reionization occurs a t a redshift near the upper end of the range 
suggested by the WMAP polarization results (ISpergel et al.l 120071 ). To explore how the 
presence of a UV background will influence our conclusions, we have run simulations with 
UV field strengths of J21 = 0.01 and 0.1, as well as with no background present. 

Simulations were run for ~ 1 Hubble time, tn, which at 2; = 25 corresponds to about 
100 Myr. At times t > tn, individual dark matter halos can no longer be treated as isolated 
objects, since most halos will have undergone at least one major merger or interaction by 
this time. For our fiducial halo, the free-fall time tg ~ (Gp)^^^^ varies between 4 Myr near 
the central cusp and 450 Myr in the outer envelope. 

For ease of reference, we have given the runs discussed in the following sections simple 
alphanumeric identifiers. We have (arbitrarily) divided up our model halos into three size 
classes: "small" halos, with masses M < 10^ Mq; "intermediate" halos, with 10^ < M < 
4.5 X 10^ Mq] and "large" halos with M > 4.5 x 10^ Mq. The first character of the identifier 
indicates with which of these classes the run in question is associated: "S" denotes a small 
halo, "I" an intermediate halo, and "L" a large halo. The next two digits in the identifier 
denote the initial redshift of the simulation. Next, there follows a character denoting the 
metallicity: "Z" corresponds to zero metallicity, while "L," "M," "N," and "S" correspond to 
metallicities 10~^, 10~^, 10~^, and 1 Zq, respectively. Following this is a number denoting the 
strength of the UV background: 1 corresponds to no background, while 2 and 3 correspond 
to J21 = 0.01 and 0.1, respectively. Finally, runs in which dust is included have a "D" 
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appended. To give an example: run I25-M3D is a run performed at redshift z = 25 with an 
intermediate-sized halo that has a metallicity Z = 10~^ Zq, an ultraviolet background field 
strength of J21 = 0.1, and that includes the effects of dust. 



Results 



3.1. Zero-Metallicity Gas 

Before examining the behaviour of the metal-enriched gas, it is prudent to first examine 
what happens in the absence of metals, i.e. in gas of wholly primordial composition. Previous 
work suggests that there is a critical protogalactic mass Merit such that gas in protogalaxies 
with M > Merit is able to cool an d collapse within a Hubble time, while gas in protogalaxies 



with M < Merit cannot (see e.g. iTegmark et al.l 119971: iMachacek et al.ll200ll ; lYoshida et al. 
2OO3I : iMashchenko et allbood : lO'Shea fc Norman! boOTh . 



To verify that we recover the same behaviour, we performed a set of three runs at 2; = 30 
with halo masses spanning a range of 2 orders of magnitude. The halos in these three runs, 
which we designate as S30-Z1, I30-Z1, and L30-Z1, had masses of 4.0 x 10^ M©, 4.8 x 10^ Mq 
and 4.8 x 10^ Mq respectively. In Figure [21 we show the time evolution of the mean number 
density and mean temperature within the scale radius in these runs. 

Table [6] gives the properties of the gas within the scale radius after 100 Myr of evo- 
lution. The table lists the minimum central temperature Tc^min and the maximum central 
density nc,max- We also calculate the fractional H2 abundance with respect to the total 



number of hydrogen nuclei, and list the maximum central value, Xh 



2,c,max; 



in the table. 



The mass fraction of gas within a sink particle is denoted by /sink- Because we have no 
information on the density and temperature distribution within the sink particle radius, we 
adopt the values at time of sink particle formation in Table [61 Note that as the gas in run 
S30-Z1 does not collapse, it never reaches the density required for a sink particle to form. 

Both Figure [H and Table M clearly demonstrate the change in behaviour that occurs for 
increasing M: gas in run L30-Z1 cools and collapses very quickly, on a timescale t <^ t^, while 
gas in run I30-Z1 cools and collapses on a slower timescale t ~ tn- Finally, the gas in run 
S30-Z1 cools, but not by enough to bring its temperature below the halo virial temperature, 
and so it does not collapse within a Hubble time. Indeed, we have run this simulation for 
~ 520 Myr longer than is shown in the plot, but find no evidence for collapse even after this 
extended period. We can therefore conclude that 4.0 x 10"^ Mq < Merit < 4.8 x 10^ Mp^, with 



the value likely somewhere near the upper end of this range. For comparison, iTegmark et al. 
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fll997h find that Mcnt = 10^ M© at 2 = 30; Machacek et al.l ( 20011) find a s i gnific antly 
smaller value of Merit = 1-25 x 10^ Mq at the same redshift; andlYoshida et al.l (120031 ) find 
that Merit ^ 7.0 X 10^ M© at 2 < 25. More recently, lO'Shea fc Norman! ((2003) haye shown 



that in their ensemble of simulations of primordial star formation, the mass of the first star- 
forming halo lies in the range 1.4 X 10^ < M < 6.9 X 10^ Mq, where the scatter appears to 
be a result of the sensitivity of the problem to the initial conditions. Our own estimate of 
Merit is in general agreement with these determinations, suggesting that our basic treatment 
of the physics of the primordial gas is sound. 

Another area that we can explore with our zero-metallicity runs is the response of the 
gas to the presence of an ultraviolet background. Previous work has shown that H2 cooling 
of protogalactic gas is significantly affected once the strength of the background exceeds 
J21 — 0.01 and that collapse in small halos is almost entirely suppressed for J21 > 0.1, 



owing to efficient photodissociation of the H 2 (iHaiman et al.l l2000l : iMachacek et al. 



up 

M 



2001 



Yoshida et al.l l2003l : iMashchenko et al.l |2006| ) . We show in Figure [3] an example of what 
happens in our own simulations as the strength of the background is increased. In the figure, 
we plot results from three runs performed at 2; = 25 with halo mass M = 7.8 x 10^ M©. These 
runs are designated I25-Z1, I25-Z2, and I25-Z3 and had UV background field strengths of 
J21 = 0.0, 0.01 and 0.1, respectively. We see from the figure that the imposed UV background 
of strength J21 = 0.01 in run I25-Z2 delays collapse by more than 30 Myr compared to run 
I25-Z1, but that collapse does eventually occur. On the other hand, a further increase in 
the UV background field strength by a factor of 10 completely inhibits the collapse of the 
gas within the time frame of our simulation. These results are in good agreement with the 
previous studies cited above, further increasing our confidence that our basic modelhng is 
sound. 



3.2. Metal-Enriched Gas 



3.2.1. Low-Metallicity Protogalaxies 



Having established how primordial gas behaves in these conditions, we next examined 
the effect of enriching the gas with metals. As previously mentioned in § 12. 3^ we focussed 
initially on gas with a metallicity Z = 10~^ Zq. This metallicity is at the upper limit of the 
range of values proposed for the so-called critical metallicity Zcrit, the value of the metallicity 
at which efficient fragmentation and low-mass star formation is hypothesized to first occur 
JSromm et al.lboOll : liromm fc Loebll2003l : ISchneider et allbooi lOmukai et al.ll2005[ l. It is 
also comparable to t he globally averaged meta llicity produced by the sources responsible for 



reionization (see e.g. iRicotti &: Ostrikerl 12004 ). 
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In the top panel of Figure H] we show the time evolution of the relative difference between 
the central density in run I25-L1 and in I25-Z1. In other words, we plot |nz c ~ '^l,c|/'^z,c; 
where nL,c and riz ^ are the central densities in runs I25-L1 and I25-Z1, respectively. We plot 
this quantity to highlight the difference between the two runs, which would be difficult to 
discern in a more conventional comparison plot. We see that prior to sink particle formation, 
the central density of the two runs differs by less than 10%, and that at early times (t < 
50 Myr), the difference is less than 1%. In the bottom panel of Figure H] we show a similar 
plot for the central temperature in the two runs. Again, the difference between the two runs 
is very small. 

These results demonstrate that the cooling of the gas at low densities (n < 1 cm~'^) and 
high temperatures (T > 2000 K) is barely influenced by the presence of metals: fine-structure 
cooling contributes only marginally to the total cooling rate, and the results of the runs agree 
to within 1%. The evolution of cooling during the first 50 Myr and the onset of collapse are 
thus almost independent of the metallicity of the gas, at least for gas with Z < 10~^ Zq. 
At later times during the collapse, the difference between the metal-enriched and the zero- 
metallicity models becomes more pronounced, indicating that metal-line cooling is beginning 
to play a more significant role. Despite this, the behaviour of the gas in the two runs remains 
very similar up to the point at which sink formation occurs (i.e. at n < 500 cm~^). 

An obvious question to ask is whether these results are somehow peculiar to the partic- 
ular combination of halo mass and redshift that we have chosen to examine. To investigate 
this, we performed additional simulations with both Z = and 10~^Zq for a range of different 
halo masses and redshifts, as summarized in Table HI We found that in all cases we obtained 
similar results; the thermal and dynamical evolution of gas in runs with Z = 10~^ Zq differed 
very little from that of gas in runs with Z = 0. 



3.2.2. The effects of rotation 



We have also investigated the possibility that our results were significantly affected by 
our neglect of the effects of rotation in the majority of our runs. Theoretically, we do not 
expect rotation to b e important o n scales r ^ Arvir, where A is the dimensionless spin 
parameter, given by (|Peebleslll97ll ) 

where J is the total angular momentum of the halo, M is the halo mass and E is the total 
(kinetic plus potential) energy of the halo. Previous work has shown that for halos with 
the range of masses and redshifts considered in this paper, A has a lognormal distribution. 
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with mean A = 0.035 flYoshida et al.ll2003l ). Therefore, we expect that rotation will only 
significantly affect the dynamics of the gas within the innermost few percent of the halo. As 
this is typically comparable to the radius at which we cease to be able to resolve the collapse 
of the gas, we do not expect the inclusion of rotation to significantly affect our results or to 
materially alter our conclusions. 

In order to verify that this is a reasonable expectation, we performed one simulation 
in which the initial gas distribution was given a non-zero angular momentum. Within the 
virial radius of the halo, the gas was placed into rotation with constant angular velocity. At 
larger radii, the initial angular velocity decreased linearly with radius, reaching zero at the 
truncation radius of the halo. This run, which we designate hereafter as I25-Z1-ROT, had 
a spin parameter A = 0.05. The other parameters used for this run were the same as in 
run I25-Z1. In Figure we compare the outcome of these two simulations. We see that, as 
expected, rotation has very little effect on the evolution of the gas at early times. It does 
appear to affect the evolution once the density exceeds n = 100 cm~^, significantly slowing 
the collapse, but at this point we are close to the resolution limit of our simulation and so 
we are unable to say with any certainty whether or not a rotationally supported disk forms. 
Nevertheless, as the main results of this paper concern the evolution of the gas at densities 
n -C 100 cm~^, it is clear that they are unaffected by our neglect of rotation in the majority 
of our simulations. 



3.2.3. When Do Metals Make a Difference? 

If cooling from metals is ineffective when Z = 10^^ Zq, then just when does it be- 
come important? How enriched does the gas need to be before metal-line cooling makes a 
significant difference? To investigate this, we performed a set of additional runs using our 
fiducial halo parameters {z = 25, M = 7.8 x 10^ Mq) that looked at the effects of further 
increasing Z. In Figure O we plot the time evolution of the central density and central 
temperature in runs I25-Z1, I25-L1, I25-M1, I25-N1, and I25-S1, which had metallicities 
of Z = 0.0, 10"'^, 10~^, 0.1 and 1.0 Zq, respectively. We also show in the plot the results 
of runs I25-L1D, I25-M1D, and I25-S1D, which had metallicities of Z = 10^^^ 10"^ and 
1.0 Zq, respectively, and which differed from runs I25-L1, I25-M1, and I25-S1 by including 
the chemical and thermal effects of dust. 

It is clear from the figure that all of the runs with Z < 0.01 Zq evolve in essentially 
the same way. Some minor differences in the rate of collapse and in the central temperature 
of the gas are apparent at late times, particularly in the two runs with Z = 0.01 Zq, but 
these do not appear likely to significantly affect the outcome of the simulation. If we enrich 
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the gas by another order of magnitude, to Z = 0.1 Zq, however, then we begin to see the 
metals playing a significant role in the cooling of the gas. In run I25-N1 the gas cools more 
rapidly than in any of the lower metallicity runs, and as a consequence, runaway gravitational 
collapse sets in approximately lOMyr earlier than in the lower metallicity runs. Nevertheless, 
the evolution of the gas in this run remains qualitatively very similar to that in the lower 
metallicity runs. It is only when we enrich the gas by a further order of magnitude, to solar 
metallicity, that we see a truly dramatic difference in the outcome of the simulation, with 
the gas in runs I25-S1 and I25-S1D cooling and collapsing within the first 20Myr of the run. 
Based on these results, we would estimate the critical metallicity at which metal cooling 
dominates in the low-density gas to be Z^rit ~ 0.1 Zq, to the nearest order of magnitude. 

Figure M also allows us to assess the role played by dust in these systems. The fact that 
runs I25-L1D and I25-M1D produce essentially identical results to runs I25-L1 and I25-M1 
demonstrates that at Z = 0.01 Zq and below, dust has a minimal effect on the low-density 
evolution of these systems. On the other hand, at high metallicities, dust does significantly 
affect the evolution of the gas. Although the presence of dust might be expected to aid the 
cooling of the gas, it is clear from Figure [6] that in actual fact the gas cools more rapidly 
when dust is not present. This appears to be a consequence of the fact that the ionized gas 
recombines more rapidly when dust is present, owing to recombination of ionized hydrogen 
on the surface of dust grains, with the result that both the Compton cooling rate, as well 
as the contributions of electron-ion and electron-atom collisions to the fine-structure cooling 
rate, fall off more rapidly with time. 



3.2.4- Collapse with a UV background 

Finally, it is important to determine whether metal-line cooling allows collapse to occur 
in UV-irradiated halos that would otherwise be unable to collapse. Figure [7] summarizes the 
results of a set of simulations designed to address this question. It shows the evolution of 
the central density and temperature of the gas in a set of runs in which J21 = 0.1, M and 
z were held constant at their fiducial values, and the metallicity was varied. Details of the 
eventual state of the gas at the end of these runs are given in Table [61 We see that even at 
metallicities as high as Z = 0.01 Zq, an ultraviolet background flux of strength J21 = 0.1 
is sufficient to completely suppress collapse. We therefore conclude that in low-metallicity, 
low-mass protogalactic halos, metal-line cooling cannot initiate collapse. Either the surviving 
H2 provides enough cooling to allow collapse to occur, irrespective of whether or not there 
are metals present, or the gas does not collapse. Only at quite high metallicities {Z > 0.1 Zq) 
do the metals actually make an appreciable difference to the evolution of the gas. 
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4. Discussion 

At first sight, tlie fact tliat fine structure cooling from metals has little impact on 
the thermal or dyna mical evolution of the gas at metallicities below 0.1 Zq is somewhat 



surprising, given that iBromm et al.l (120011 ) found that gas with a metallicity of only Z = 
1O~'^Z0 could cool rapidly and fragment even in the complete absence of molecular hydrogen. 
However, comparison of the cooling time due to fine-structure emission with the free-fall time 
of the gas helps to make the situation clear. The cooling time due to fine- structure emission 
is given by 

X Tl/Jv'J- ' 

''^cool.fs = T—: , (3) 

7-1 Afs 

where Afg is the total cooling rate per unit volume due to fine-structure emission, T is the 
temperature of the gas, 7 is the adiabatic index, n is the number density of the gas, and k 
is the Boltzmann constant. The free-fall time can be written as 

where G is the gravitational constant and p = Pgas + Pdm- The gas density pgas can be written 
in terms of the number density n as Pgas = pn, where is the mean molecular weight of the 
gas. 

In Fi gure [HI we indicate the temperatures and densities at which tcooi,fs — '^ff for 10 Zq 
gas, for three different assumed fractional ionizations, with Xe = 1.0 (fully ionized gas), 10^^, 
and 10~^. In each case, we assume that the carbon and silicon are present only as C+ or Si^, 
since a fairly small external UV flux is sufficient to achieve this. In the Xe = 1.0 case, we also 
assume that all of the oxygen is O"*", since charge transfer between oxygen and hydrogen, 
which have nearly identical ionization potentials, ensures that xo+/xo — xu+fx-a. For the 
dark matter density pdm we adopt the value found at the center of our fiducial protogalactic 
halo at t=0, pdm — 8 x 10~^^gcm~^. In the plot, regions to the left of the line have tcooi,fs > ts, 
while those to the right have tcooi,fs < tg. 

At the beginning of our fiducial simulation, gas in the center of the halo has a temper- 
ature T = 10^ K and a number density n = 0.03 cm~^. It therefore lies outside of the regime 
where fine-structure cooling is efficient, and so it is not surprising that we find that metal-line 
cooling is initially unimportant. As the gas cools, whether through Compton cooling or H2 
emission, and begins to compress as it falls into the halo, it moves towards the temperature 
and density regime in which fine-structure cooling is effective. As an example, we plot in 
Figure [8] the trajectory in the density-temperature plane followed by gas at the center of the 
halo in run I25-L1. However, at the same time as the gas is cooling and collapsing, it is 
also recombining. This causes the boundary of the efficient cooling regime to move to the 
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right of the plot, towards higher densities. The physical reason for this shift is the fact that 
free electrons are much more effective than neutral hydrogen at exciting the C"*" and Si"*" 
fine-structure lines, and so electron excitation dominates for Xe > 10~^-10~'^, depending on 
the temperature. The net effect is that fine-structure cooling remains of little importance 
until the gas is near the high-cooling regime. This only happens after considerable cooling 
and compression has already taken place, and therefore does not occur at all if H2 cooling is 
ineffective, as is the case, for instance, in runs I25-Z3 or I25-L3. 

For gas with a metallicity greater than 10~^ Zq, fine-structure cooling is effective in a 
larger region of the density-temperature plane. At low densities [n < 1 cm~^), an increase 
in Z of an order of magnitude shifts the solid lines in Figure [8] by an order of magnitude 
with respect to density. This allows us to estimate the metallicity at which fine-structure 
cooling will dominate even in the initial low-density gas. We find that fine-structure cooling 
should dominate for metallicities in excess of 0.1-1.0 Zq, depending on the details of the 
evolutionary trajectory followed by the gas at early times. This estimate agrees well with 
the results we obtain from our SPH simulations. 

Changes to the halo mass or redshift of collapse will affect this estimate in two ways. 
First, increasing the redshift of collapse will make the dark matter halo denser and will 
therefore decrease tg at early times, when dark matter is gravitationally dominant. This will 
shrink the region in which fine-structure cooling can effectively operate. Second, altering the 
halo mass and redshift will alter the evolutionary trajectory followed by the gas. Specifically, 
since the initial temperature of the gas is fixed at 10^ K in these simulations, changes in M 
or z that increase the virial temperature of the halo will make the gas more gravitationally 
unstable, allowing it to collapse more quickly. The collapsing gas will therefore be hotter 
than in our fiducial case (as it will have had less time to cool), and will also be slightly more 
ionized. On the other hand, if the virial temperature of the halo is lowered, then the gas will 
be more stable, will collapse more slowly, and will therefore be cooler and less ionized than 
in our fiducial case. However, in either case, we would not expect to see a change of more 
than a factor of a few in the critical metallicity required for effective fine-structure cooling 
at low densities, and our estimate should therefore still be of the correct order of magnitude. 

Given these results, we next ask whether they can possibly be consistent with the results 



of iBromm et al.l (120011 ). who find that in their simulations, fine structure cooling is effective 
even for Z = 10~^ Zq, which is at least 2 orders of magnitude smaller than our value. In fact, 
their results do actually appear to be consistent with ours, despite the apparent disagreement. 
The reason is due to the difference in the initial conditions used for their simulations and 
for ours. The gas in their simulations has an initial temperature of 200 K aX z = 100, and 
its temperature falls further due to adiabatic cooling in the IGM prior to the formation of 
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their simulated protogalactic halo at 2; ~ 30. Since this halo has a mass of 2 x 10^ Mq and 
a virial temperature T^i-^ ~ 5000 K, the temperature of the gas at the moment at which the 
halo forms is very much smaller than the halo virial temperature. Consequently, thermal 
pressure is initially unable to prevent the collapse of gas into the halo. As it collapses, the 
gas heats up, and the thermal pressure eventually becomes large enough to halt the collapse. 
However, this does not occur until the gas temperature is approximately equal to the halo 
virial temperature, by which time the gas density has increased to a value n ~ 300 cm~^. 
As can be seen from Figure [SI gas with this density and with a temperature of 5000 K lies 
close to or within the regin ae where tonoi.fa < tff (depending on its fractional ionization), and 
so it is not surprising that iBromm et al.l (120011 ) find that fine-structure cooling is effective 
and that the gas ca r i cool even in the complete absence of H2. It is also important to stress 
that iBromm et al.l (I2OOII ) do not find fine-structure cooling to be important at densities 
n -C 100 cm~^, as the almost adiabatic initial evolution of the gas makes clear. 

For comparison, our simulations start with a high initial gas temperature of 10^ K and 
involve halos with virial temperatures Tvir < 10"^ K; for instance, Tvir — 1900 K in our fiducial 
simulation. This means that in our simulations, Tg^s > T^w initially, with the result that 
thermal pressure support is important right from the start. Indeed, at the beginning of 
our simulations pressure precisely bala nces gravity, by desi gn. Therefore, there is no initial 
phase of pure free-fall collapse as in the lBromm et al.l (I2OOII ) simulations. Instead, significant 
gravitational contraction occurs only if the gas is able to cool to a temperature of order Tvir 
or below, which, since fine-structure cooling is initially ineffective, will only occur if enough 
H2 can form in the low-density gas. At late times in our simulations, we wo uld expect to 



see fin e-structure cooling become more effective at lower Z, in line with what iBromm et al. 



( I2OOII ) find, but unfortunately the required gas density is only marginally resolvable with 
our current simulations. 

While the two sets of simulations therefore appear to be consistent with each other, a 
key question is obviously which set of initial conditions is more appropriate for describing 
the early evolution of small, metal enriched protogalaxies. We argue that it is difficult to 
see how intergalactic gas could become metal enriched without at some point being ionized. 
Previous calculations have shown that the size of a typical region enriched by a Population 
III super nova is much smalle r than the size of the H ll region created by its progenitor star 
(see e.g. iBromm et al.ll2003l). We would therefore expect our initial conditions to be more 



appropriate than those of iBromm et al.l (120011 ) for treating recently enriched and ionized 
regions. However, if enough time elapses following the enrichment event for the gas to be 
able to cool down to a temperature of a fe w hundred K, then th e Bromm et al. (2001) initial 
conditions will be more appropriate. As lOh fc HaimanI (120031 ) show, this is most likely to 
occur in high-redshift gas with a low overdensity. Then Compton cooling is fast and highly 
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effective and can cool the gas to T ~ 300 K within a recombination time. On the other 
hand, at lower redshifts, or at higher overdensities, the gas recombines faster than it cools, 
and the temperature that can be reached by Compton cooling alone is much higher. This is 
the case in our simulated halos. Note that in either case the metals play no significant role 
in determining whether or not the gas is able to collapse. 

An important implication of these results is that if we are primarily concerned with 
investigating questions such as how M^rit evolves with redshift, or how UV feedback in the 
form of Lyman- Werner photons affects the ability of the gas to cool, then we need not worry 
about the effects of metal enrichment. This is because the thermal evolution of the gas on the 
scales of interest for these questions is complet ely dominated by Com pt on cooling and /or Ho 



cooling. Therefore, results from studies such as lHaiman et al.l (120001 ) or lYoshida et al.l (120031 ) 



give a better guide to the behaviour of small, low-metallicity protogalaxies than might have 
been anticipated (although the additional complications posed by the mechanical energy 
injected into the gas by H II regions and supernovae do of course still need to be taken into 
account). 

Finally, it is important to stress that our results do not address the question of whether 
or not there is a critical metallicity Zait above which fine structure cooling from metals 
allows efficient fragmentation to occur. This is because if fragmentation does occur, we 
would expect it to occur at densities n > 500 cm~^, which are unresolved in our current 
simulations. We will examine this question with much higher resolution simulations in a 
future paper. 
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Table 1. List of Reactions Included in Our Chemical Network 

Number Reaction Reference 

1 H + e" -> H- + 7 Wishart (1979) 

2 H- + H ^ H2 + e" Launav et al. (1991) 

3 H + H+ ^ H+ + 7 Ramaker fc Peek (1976) 

4 H + H+ ^ H2 + H+ Karpas et al. (1979) 

5 H- + H+ ^ H + H Moselev et al. (1970) 

6 H+ + e- ^ H + H Schneider et al. (1994) 

7 H2 + H+ ^ H+ + H Savin et al. (2004) 

8 H2 + e" H + H + e" Stibbe &: Tennyson (1999) 

9 H2 + H^H + H + H Mac Low fc ShuU (1986) 

10 H2 + H2 ^ H2 + H + H Martin et al. (1998) 

11 H + e" ^H+ + e- + e- .Tanev et al. (19871 
13 H+ + e' ^ H + 7 Ferland et al. (1992) 

15 H~ + e" H + e" + e" Janev et al. (1987^1 

16 H- + H ^ H + H + e" Janev et al. (1987) 

17 H- + H+ ^ H+ + e- Poulaert et al. (1978) 
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Table 1 — Continued 

Number Reaction Reference 

30 C+ + e- ^ C + 7 Nahar fc Pradhan (1997} 

31 Si+ + e" ^ Si + 7 Nahar (20001 

32 0+ + e" ^ O + 7 Nahar f 19991 

33 C + e- ^ C+ + e- + e- Voronov (1997') 

34 Si + e" ^ Si+ + e" + e' Voronov (1997) 

35 O + e" ^ 0+ + e- + e" Voronov (1997) 

36 0+ + H ^ O + H+ Stancil et al. (1999) 

37 O + H+ ^ 0+ + H Stancil et al. (1999) 
39 C + H+ ^ C+ + H Stancil et al. (1998) 

42 Si + H+ Si+ + H Kingdon fc Ferland (1996) 

44 C+ + Si ^ C + Si+ Le Teuff et al. (20001 

51 H- + 7 ^ H + e" Wishart (1979) 

52 H+ + 7 ^ H + H+ Dunn (1968) 

53 H2 + 7 H + H Draine fc Bertoldi (1996) 
56 C + 7 C+ + e" Verner et al. (1996} 

58 Si + 7 Si+ + e" Verner et al. f 19961 
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Table 1 — Continued 



Number 


Reaction 




Reference 


60 


H + H(+grain) — 




HoUenbach & McKee fl979) 


61 


H"*" + e~(+grain) 


^ H 


Weineartner & Draine f2001') 


64 


C"*" + e~(+grain) 


^ C 


Weinecartner & Draine (2001) 


66 


Si"*" + e~(+grain) 


Si 


Weineartner & Draine f 20011 



Note. — The reaction numbering scheme used here is the same 
as that in paper I. References are to the primary source of data for 
each reaction. 
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Table 2. List of Photochemical Reaction Rates. 



Number 


Reaction 




Rate {J 21 s"^) 




51 


H-+7- 


H + e- 


i?5i = 1-36 X 10- 


ll 


52 


H+ + 7- 


H + H+ 


i?52 = 4.11 X 10- 


-12 


53 


H2 + 7- 


> H + H 


i?53 = 1.30 X 10- 


-12 


56 


C + 7^ 


C+ +e- 


i?56 = 5.56 X 10- 


-12 


58 


Si + 7-. 


Si+ + e~ 


i?58 = 2.44 X lo- 


-11 



Note. — The reaction numbering scheme used 
here is the same as that in Paper I. Rates are cal- 
culated assuming an incident spectrum correspond- 
ing to a modified, diluted 10^ K black body, as de- 
scribed in the text. The value J21 quantifies the 
strength of the radiation field at the Lyman limit: 
J{i-'a) = 10^^^ J21 erg s~^ cm~^ Hz^^ sr^^. The rate 
listed for H2 photodissociation (reaction 53) is for 
unshielded gas; our treatment of H2 self-shielding is 
described in section 12. 2[ 
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Table 3. Processes Included in Our Thermal Model. 



Process 



References 



Cooling 

Fine-structure lines (C, C+, O, Si, Si^ 
Lja cooling 
H coUisional ionization 
Compton cooling 
recombination 
H2 rovibrational lines 
H2 collisional dissociation 
Gas-grain energy transfer" 
Heating 

Photoelectric effect 

H2 photodissociation 

UV pumping of H2 

H2 formation on dust grains 



See Paper I 
Cen (1992) 
Janev et al. (1987) 
Cen (1992) 
Ferland et al. (1992) 
Le Bourlot et al. (1999) 

Mac Low fc Shull f 19861 : Martin et al. fl998) 
HoUenbach fc McKee (1989) 

Bakes fc Tielens fl994l : Wolfire et al. fl995l 
Black fc Dakarno (1977) 
Burton et al. (1990) 
HoUenbach fc McKee (1989) 



Note. — a: If T„a<^ < T, 



grain ; 



the net flow of energy is from the grains to the gas, leading 



to heating instead of cooling. 
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Table 4. List of the Runs Discussed in this Paper 



Run 


Mass(M0) 


Redshift 




Dust^ 




S30-Z1 


4.0 


X 


10* 


30 


0.0 


No 


0.0 


I30-Z1 


4.8 


X 


10^ 


30 


0.0 


No 


0.0 


I30-L1 


4.8 


X 


10^ 


30 


10-3 


No 


0.0 


I25-Z1-R0T^ 


7.8 


X 


10^ 


25 


0.0 


No 


0.0 


I25-Z1 


7.8 


X 


10^ 


25 


0.0 


No 


0.0 


I25-L1 


7.8 


X 


10^ 


25 


10-=^ 


No 


0.0 


I25-L1D 


7.8 


X 


10^ 


25 


10-=^ 


Yes 


0.0 


I25-M1 


7.8 


X 


10^ 


25 


10-2 


No 


0.0 


I25-M1D 


7.8 


X 


10^ 


25 


10-2 


Yes 


0.0 


I25-N1 


7.8 


X 


10^ 


25 


10-1 


No 


0.0 


I25-S1 


7.8 


X 


10^ 


25 


1.0 


No 


0.0 


I25-S1D 


7.8 


X 


10^ 


25 


1.0 


Yes 


0.0 


I25-Z2 


7.8 


X 


lO'^ 


25 


0.0 


No 


10-2 


I25-Z3 


7.8 


X 


10^ 


25 


0.0 


No 


10-1 


I25-L3 


7.8 


X 


10^ 


25 


10-=^ 


No 


10-1 


iZO-iViO 


7.8 


X 


10^ 


zo 


lU 


IMO 


lU 


I25-M3D 


7.8 


X 


10^ 


25 


10-2 


Yes 


10-1 


I25-N3 


7.8 


X 


10^ 


25 


10-1 


No 


10-1 


I25-S3 


7.8 


X 


10^ 


25 


1.0 


No 


10-1 


I20-Z1 


1.5 


X 


10*^ 


20 


0.0 


No 


0.0 


I20-L1 


1.5 


X 


10^ 


20 


10-3 


No 


0.0 


L30-Z1 


4.8 


X 


106 


30 


0.0 


No 


0.0 


L30-L1 


4.8 


X 


10^ 


30 


10-3 


No 


0.0 


L20-Z1 


1.2 


X 


10^ 


20 


0.0 


No 


0.0 


L20-L1 


1.2 


X 


10^ 


20 


10-3 


No 


0.0 



^Indicates whether the run includes dust or only gas-phase 
metals 

^Strength of the UV background 
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^This run included initial rotation, with a spin parameter 
A = 0.05 
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Table 5. Halo Mass Resolution and Virial Temperature 



Halo mass (Mq) 


Redshift 


Particle number N 








4.0 X 10^ 


30 


132651 


0.25 


20 


310 


4.8 X 10^ 


30 


117649 


2.0 


160 


1600 


7.8 X 10^ 


25 


140608 


2.5 


200 


1900 


1.5 X 10^ 


20 


140608 


5.0 


400 


2400 


4.8 X 10*^ 


30 


493039 


5.0 


400 


7600 


1.2 X 10^ 


20 


1042189 


5.0 


400 


9500 
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Table 6. Physical State of the Densest Gas within the Scale Radius at Time tend 



Run 




T b 
-'21 


^end 


rp d 
c,min 


"'c,max 


•^H2,c,max 








(Ze) 




(Myr) 


(K) 


(cm-3) 








S30-Z1 


0.0 


0.0 


100 


1100 


0.3 


9 X 10-^ 


0.0 




I30-Z1 


0.0 


0.0 


75 


< 200 


> 500 


5 X 10-3 


0.04 




I30-L1 


10-3 


0.0 


100 


< 200 


> 500 


3 X 10-3 


0.19 




I25-Z1-ROT 


0.0 


0.0 


100 


< 200 


> 500 


4 X 10-3 


2 X 10- 


-3 


I25-Z1 


0.0 


0.0 


100 


< 200 


> 500 


4 X 10-3 


0.17 




I25-L1 


10-3 


0.0 


100 


< 200 


> 500 


2 X 10-3 


0.20 




I25-L1D 


10-3 


0.0 


100 


< 200 


> 500 


0.5 


0.13 




I25-M1 


10-2 


0.0 


100 


< 200 


> 500 


3 X 10-3 


0.15 




I25-M1D 


10-2 


0.0 


100 


< 200 


> 500 


0.5 


0.08 




I25-N1 


10-1 


0.0 


100 


< 200 


> 500 


2 X 10-3 


0.17 




I25-S1 


1.0 


0.0 


100 


< 200 


> 500 


1 X 10-3 


0.28 




I25-S1D 


1.0 


0.0 


100 


< 200 


> 500 


0.5 


0.27 




I25-Z2 


0.0 


10-2 


100 


< 200 


> 500 


5 X 10-^ 


5 X 10" 


■3 


I25-Z3 


0.0 


10-1 


100 


4400 


0.7 


1.6 X 10-5 


0.0 




I25-L3 


10-3 


10-1 


100 


5900 


0.6 


1.5 X 10-5 


0.0 






in-2 

lU 




1 nn 

lUU 




n R 

U.O 




n n 




I25-M3D 


10-2 


10-1 


100 


5900 


0.7 


1.5 X 10-5 


0.0 




I25-N3 


10-1 


10-1 


100 


2600 


2.5 


5 X 10-5 


0.0 




I25-S3 


1.0 


10-1 


100 


< 200 


> 500 


5 X 10-^ 


0.3 




I20-Z1 


0.0 


0.0 


100 


160 


1600 


3 X 10-3 


0.0 




I20-L1 


10-3 


0.0 


100 


130 


2100 


3 X 10-3 


0.0 




L30-Z1 


0.0 


0.0 


25 


< 200 


> 500 


7 X 10-3 


1 X 10" 


-3 


L30-L1 


10-3 


0.0 


10 


< 200 


> 500 


3 X 10-3 


2 X 10- 


-4 


L20-Z1 


0.0 


0.0 


20 


< 200 


> 500 


7 X 10-3 


6 X 10- 


-4 


L20-L1 


10-3 


0.0 


10 


310 


330 


3 X 10-3 


0.0 





*MetaUicity of the gas. 
''Strength of the UV background. 



-32- 



'^Time at the end of the simulation, or t = 100 Myr, whichever is smaller 
^Minimum temperature of the gas within the scale radius r^. 
^Maximum number density of the gas within the scale radius rg- 
^Maximum fractional H2 abundance within the scale radius r^. 
^Mass fraction of gas within a sink particle. 



-33- 




Temperature (K) 

Fig. 1. — Comparison of various parameterizations of the H2 cooling function, plotted in units 
of ergcm^s"^. Rates are computed for a number density n = Icm^'^ and assuming that uh 3> 
and that the H2 ortho-to-para ratio is 3:1. Solid line: Le Bourlot et al. (1999); dashed 
line: Galli & Palla (1998); dotted line: Lepp & ShuU (1983); dash-dotted line: HoUenbach & 
McKee (1979). 
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100 - 




100 I 1 1 1 1 1 

20 40 60 80 100 

Time [Myr] 

Fig. 2. — Top: Time evolution of the gas density within the scale radius of the dark matter 
halo in runs L30-Z1 (solid line), I30-Z1 [dashed line), and S30-Z1 (dotted line), corresponding 
to halos with masses 4.8 X 10^ 4.8 X 10^ and 4.0 x 10"^ Mq, respectively. Bottom: bame as 
the top panel, but for the central temperature of the gas. 
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20 40 60 80 100 

Time [Myr] 

Fig. 3. — Top: Same as the top panel of Fig. [2], but for runs I25-Z1 {solid line), I25-Z2 
{dashed line), and I25-Z3 {dotted line). These runs were started at z = 25 with a halo 
mass M = 7.8 X 10^ Mq, and had UV background field strengths J21 = 0.0, 0.01, and 0.1, 
respectively. Bottom: Same as the top panel, but for the central temperature of the gas. 
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20 40 60 80 

Time [Myr] 



Fig. 4. — Top: Time evolution of the relative difference \nz,c — ^t-LjcI/^z.c, where nL,c is the 
central gas density in run I25-L1 and nz,c the central density in run I25-Z1. For clarity, wc 
only plot the evolution until the point at which a sink particle forms (or until the end of 
the run, if no sink forms). Bottom: Same as the top panel, but for the relative difference 
between the central temperature in the two runs. 
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20 40 60 80 100 

Time [Myr] 

Fig. 5. — Top: Time evolution of the gas density within the scale radius Vg of the dark 
matter halo in runs I25-Z1 {solid line) and I25-Z1-ROT {dashed line). Bottom: Same as the 
top panel, but for the central temperature of the gas 
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10 20 30 40 50 60 70 80 

Time [Myr] 

Fig. 6. — Top: Time evolution of the gas density within the scale radius of the dark matter 
halo in runs I25-Z1 {solid line), I25-L1 {dot-dot-dashed line), I25-M1 {dash-dotted line), 125- 
Nl {dashed line), and I25-S1 {dotted line), which have metallicities Z — 0.0, 10~^, 10~^, 0.1, 
and 1.0 Zq, respectively. We also include three runs with dust, I25-L1D, I25-M1D, and 125- 

SID, denoted by the lines with symbols. In all of these runs, the halo mass M — 7.8 x 10^ Mq 
and the run begins at a redshift z = 25. No ultraviolet background is present in any of these 
runs. Bottom: Same as the top panel, but for the central temperature of the gas. 
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Fig. 7. — Top: Same as Fig. [6], but for runs I25-Z3 {solid line), I25-L3 {dot- dot- dashed 
line), I25-M3 {dash-dotted line), I25-N3 {dashed line), and I25-S3 {dotted line). The UV 
background field strength in all of these runs was J21 = 0.1. Bottom: Same as the top panel, 
but for the central temperature of the gas. 



*cool, fs > *ff 



101 



10-2 10-1 10° 10^ 102 103 

Number Density [cm~^] 



Fig. 8. — Solid contours: Temperature and density at which the coohng time due to fine- 
structure emission, tcooi,fs) equals the free-fall time, tg, for gas with a metallicity Z — 10~^ Zq 
and with fractional ionizations Xe — 1.0, 10~^, and 10"^. To the left of each fine, icooi.fs > ^fr, 
so metal cooling is inefficient. In every case, we assume that all of the carbon and silicon is 
present as C"*" and Si"*", respectively. In the Xe = 1.0 case, we assume that all of the oxygen 
is present in the form of O"*", but otherwise that it is all O. To compute the free-fall time, 
wc take the density to be the sum of the gas density pg and the dark matter density at the 
center of the halo pdm- The figure also shows how the temperature and density of the gas 
at the centre of the halo evolve in runs I25-L1 {solid line), I25-L2 {dashed line), and I25-L3 
{dotted line). The evolution of the fractional ionization in these runs is indicated by the star 
symbols: Xe — 1.0 for T — 10^ K, and decreases by a factor of 10 between each successive 
star. 



